The vacuum revealed: the final state of vacuum instabilities in compact stars 
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Quantum fields in compact stars can be amplified due to a semiclassical instability. This generic 
feature of scalar fields coupled to curvature may affect the birth and the equilibrium structure of 
relativistic stars. We point out that the semiclassical instability has a classical counterpart, which 
occurs exactly in the same region of the parameter space. For negative values of the coupling 
parameter the instability is equivalent to the well-known "spontaneous scalarization" effect: the 
plausible end-state of the instability is a static, asymptotically flat equilibrium configuration with 
nonzero expectation value for the quantum fields, which is compatible with experiments in the 
weak-field regime and energetically favored over stellar solutions in general relativity. For positive 
values of the coupling parameter the new configurations are energetically disfavored, and the end- 
point of the instability remains an open and interesting issue. The vacuum instability may provide a 
natural mechanism to produce spontaneous scalarization, leading to new experimental opportunities 
to probe the nature of vacuum energy via astrophysical observations of compact stars. 
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One of the great achievements of quantum field theory 
in curved spacetime is Hawking's semiclassical prediction 
of black hole evaporation. This process is made possible 
by the special nature of the vacuum in quantum field 
theory, but it is extremely feeble by any astrophysical 
standards. It was recently shown that the vacuum in 
strongly curved spacetimes might play an important role 
even in the absence of horizons [3, , so that the forma- 
tion and evolution of relativistic stars could be affected 
by semiclassical effects. 

In particular, Lima, Matsas and Vanzella (henceforth 
LMV) [l[ showed that the vacuum expectation value of 
nonminimally coupled scalar fields can grow exponen- 
tially in relativistic stars. An understanding of the final 
state of the LMV instability is important: if this final 
state supports nonvanishing scalar fields, the semiclassi- 
cal amplification of vacuum energy could have cosmolog- 
ical and astrophysical implications. 

In this paper we show that, for certain values of the 
coupling parameter, the most likely end-state of the in- 
stability is an asymptotically flat stellar solution with 
nonvanishing scalar field, which is compatible with grav- 
itational experiments in the weak-field regime. We also 
point out that these new solutions correspond to the well- 
known "spontaneous scalarization" phenomenon Q 
and that they are energetically favored over stellar so- 
lutions in general relativity (GR). Our findings support 
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the relevance of vacuum amplification scenarios. Vacuum 
amplification should be investigated carefully when de- 
vising strong-field tests of GR, for example in the context 
of gravitational collapse and gravitational- wave emission. 



I. SETUP: STATIC SOLUTIONS 

We look for static, spherically symmetric equilibrium 
solutions of the field equations admitting a nonzero scalar 
field with metric 

ds 2 = -f{r)di 2 + {\ - 2ro(r)/r) _1 dr 2 +r 2 d6 2 +r 2 sin 2 6 d(f> 2 . 

Following LMV we study a nonminimally coupled 
scalar field in the presence of a perfect-fluid star in GR. 
We consider the action 

S = J d^x^—gR + J d 4 x^—g C m , (1) 

where C m = £ sca iar + ^perfect fluid- The Lagrangian for a 
scalar field with conformal coupling £ is given by 

Scalar = -£ R ^ ~ 9^,^ ~ , (2) 

where y, denotes the scalar field mass. Here we set /i = 0; 
we will consider the massive case in a follow-up paper. 
For £ = 1/6, /i = the action is invariant under con- 
formal transformations (g^ — > Sl 2 .g M „ , $ — > S!~ 1< I>). For 
£ = 0, fj, = one recovers the usual minimally coupled 
massless scalar. The above Lagrangian corresponds to a 
viable theory of gravity, passing all weak-field tests 0, @ . 



2 




0.15 



0.10 



0.05 



0.00 




FIG. 1. Left panel: Existence diagram for constant density stars, showing the values of the ratio M/R s of uniform density 
compact objects which support a dynamical scalar field. The dash-dotted line (red in the online version) corresponds to the 
analytical prediction in the Newton ian limit, Eq. ff|). The vertical dashed line indicates the conformal-coupling value f = 1/6. 
We present this plot in "standard" geometrical units so that the compactness of our models can be easily compared to the 
maximum compactness for constant-density stars in GR: M/R s < 4/9 ~ 0.444 (Buchdahl's limit). Modulo this trivial unit 
conversion, the shaded regions in the left panel match exactly the shaded regions in LMV's Fig. 1 (see main text). In the 
upper-right region of the diagram, different critical lines correspond to "excited" equilibrium configurations and the integers 
refer to the number N of nodes in 4>(r) (see the main text). Right panel: Same for stars with a polytropic EOS. The dashed 
line (orange online) through the negative-^ region marks the configuration with maximum mass (i.e., the radial stability limit) 
for each f. 



The theory can be recast as a scalar-tensor theory of spherically symmetric solutions to the theory ([T]) with 



gravity in the Einstein frame |3| via the transformation 



In this form, both the classical counterpart of the LMV 
instability and its final state have been thoroughly stud- 
ied (Hi. 

We can explicitly construct static, asymptotically flat, 



nonzero scalar field and verify that our solutions match 
those found in 0, 0, 0], as follows. The equations of 
motion following from the Lagrangian are (hereafter we 
set c = G = 1) 



where 



(3) 
(4) 



T — 



p 2 g^ 2 - 2£$ 2 



a/3 



erfect fluid\ 



1 - 16tt^ 



(5) 



r 



We consider perfect-fluid, spherically symmetric stars 
with energy density p(r) and pressure P(r) such that 
perfect fluid = (P + p ) u " u v + g^P, where the fluid four- 
velocity = (1/y/J, 0, 0, 0). We specify some equation 
of state (EOS) P = P(p) and we impose regularity con- 
ditions at the center of the star, i.e. 

m(0) = 0, p(0) = po, $(0) = $ o , $'(0) = 0. (6) 

We also require continuity at the stellar radius R s , de- 
fined by the condition P(R S ) = 0. We focus on two dif- 



ferent stellar models: (i) the constant density stars (p = 
const) studied by LMV, and (ii) the polytropic model 
p = nm b +Kn m, b (r—l)~ 1 (n/no) r ,P = Kn mb(n/n ) r , 
with T = 2.34, K = 0.0195, m b = 1.66 x 10~ 24 g and 
no = 0.1 fm~ 3 (this is the model that was considered 
in Ref. Q in the context of spontaneous scalarization) . 
We have checked that nuclear-physics based EOS models 
would yield qualitatively similar results. 

Rclativistic stellar configurations in GR correspond to 
$ c = 0, so that $ = everywhere. For each central den- 
sity p c , we used a shooting method to search for nonzero 



3 



values of $ c such that &(r) —> as r —> oo. 

For constant-density configurations, we find solutions 
with nonzero scalar field in the shaded regions of the 
(£ , M/R s ) diagram shown in the left panel of Fig. [T] The 
right panel of Fig. Q] shows that the qualitative features 
of the existence diagram are the same for a polytropic 
EOS. These diagrams effectively reproduce previous re- 
sults obtained many years ago in the context of sponta- 
neous scalarization (see e.g. Q). It is remarkable that 
static solutions exist in the same region where the LMV 
instability operates (cf. Fig. 1 in LMV). This provides 
strong evidence that these solutions (when they are sta- 
ble) represent a plausible final state of the instability. 

In fact, the exact overlap between our own Fig. Q] and 
Fig. 1 of LMV can be proved analytically. Focus for sim- 
plicity on constant-density stars and masslcss scalar fields 
(but our reasoning applies in general). The critical lines 
in Fig. 1 of LMV represent the curves where marginally 
stable modes exist. These modes are zero-frequency so- 
lutions of Eq. (4) in LMV, where the potential is given 
by their Eq. (6). On the other hand, the critical lines 
in our Fig. [T] represent the boundary of regions where 
spherically symmetric, static solutions with nontrivial 
scalar field profiles cease to exist. These are solutions 
of the Einstein-Klcin-Gordon equations with $ = 0. As 
$ —> the spacetime becomes arbitrarily close to that of 
a constant-density star and the Klein-Gordon equation 
reduces to Eq. (4) in LMV, with potential given by their 
Eq. (6). Furthermore, the same boundary conditions ap- 
ply in both cases. Thus, the critical lines are obtained 
from the very same equations and they are indeed coin- 
cident, not just similar. 

Quite interestingly, the LMV instability threshold can 
be found analytically in the Newtonian limit M / R s <C 1 
(i.e., in the bottom left corner of Fig. [IJ. The instabil- 
ity line defines the existence of static solutions with a 
small but nonvanishing masslcss scalar field. The rele- 
vant equation in this limit is — 8ir£(p — 3-P) v E' = 0, 
where a prime denotes a derivative with respect to r, 
and we use the ansatz for the scalar field $ = (^ / 'r)e~ iut 
(i.e., we consider an s-wave). Assuming that p is constant 
in the stellar interior (this assumption holds exactly for 
uniform-density stars and is a good approximation for 
most EOSs), a regular solution at the origin and at in- 
finity that is also continuous at R s corresponds to 

24M£ = -7T 2 R s . (7) 

As shown in Fig. [TJ this prediction is in very good agree- 
ment with the LMV resultfQ. 



II. MASS AND BINDING ENERGY 

The orbits of bodies far away from the star depend 
on the star's gravitational mass M, shown as a function 
of the central baryonic density p c = m;,n(0) in the left 
panel of Fig. O In GR, maxima of this curve correspond 
to marginally stable equilibrium configurations, and all 
solutions after the first maximum are unstable to radial 
perturbations in the polytropic case (see e.g. Q). The 
baryonic mass of the configuration, 

m = nib J d 3 Xy/—gu a n(r) , (8) 

corresponds to the energy that the system would have 
if all baryons were dispersed to infinity. The normalized 
binding energy E^/M = fh/M — 1 is plotted in the inset 
of the left panel in Fig. [5] For bound (not necessarily 
stable) configurations, Eb > 0. 

A. Negative £ 

When £ < — 2, in the shaded region of Fig. Q] we found 
only a single solution coupled to a nontrivial scalar field. 
As shown in the left panel of Fig. [2] the critical central 
density for these solutions is roughly the same as in GR, 
but they have larger maximum mass than their GR coun- 
terparts (s=s 10% larger for £ < —5). Similar arguments 
to GR indicate that models after the first maximum in 
M(p c ) (lying above the dashed line in the negative-^ re- 
gion of the right panel of Fig. [1} are unstable. 

The inset of the left panel in Fig. [2] shows that similar 
deviations occur also for the binding energy. For a given 
m, the binding energy is higher than in GR, so these con- 
figurations are energetically favored over stellar solutions 
in GR. As shown in the right panel of Fig. [SJ for a fixed 
EOS the scalar field can sensibly modify the mass-radius 
relation. These modifications depend on £, but they can 
be of order 10% or more, and as such they could be 
observable. Present neutron star observations give con- 
straints on the mass-radius relationship (e.g. [10l fTl|). 
and electromagnetic observations of binaries containing 
X-ray pulsars could in principle constrain the binding en- 
ergy as we ll fl2l . [l3T | . A recent high-mass neutron star ob- 
servation ([Ij]) also rules out many EOS models in GR. If 
astrophysical measurements are not compatible with any 
realistic EOS in pure GR, we may be able to constrain 
the coupling constant £ and to probe the occurrence of 
vacuum instabilities in astrophysical environments. 



B. Positive £ 



The basic features of the instability in compact stars were under- 
stood by Ford (who studied unstable scalar fields as a possible 
mechanism to damp the effective value of the cosmological con- 
stant) as early as 1987 see also Q|. 



For £ > 1/6 the situation is more complicated. De- 
pending on the compactness M/R s and on the coupling 
constant £, we have found a hierarchy of new solutions, 
which can be labeled by the number N of radial nodes of 
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FIG. 2. Left: Gravitational mass as a function of the central baryonic density p c /po, where po = 8 x 10 14 g/cm 3 is a typical 
central density for neutron stars. The inset shows the (normalized) binding energy as a function of p c /po- Right: Gravitational 
mass as a function of the radius for different values of the coupling. 



the scalar field. In general, the larger N, the lower the 
binding energy, so high-TV solutions can be thought of as 
"excitations" of the (energetically favored) ground-state 
configuration. Critical curves bounding the shaded exis- 
tence region of these solutions are plotted in the upper- 
right corner of the two panels of Fig. [T] 

The inset of Fig. [2] shows that ground-state config- 
urations actually have a negative binding energy for 
Pc > Pcrit ~6x 10 15 g/cm 3 . When p c < p clit , we only 
found solutions with binding energy lower than in GR, 
so GR solutions are energetically favored in this region of 
parameter space. The fact that negative-^ solutions are 
more stable than positive-^ solutions is consistent with 
the interpretation of £ as an "effective gravitational con- 
stant" proposed by van der Bij and Gleiser [l5[ in their 
study of boson stars with nonminimal coupling. 

It is reasonable to conjecture that the top-right region 
of Fig. [1] does not correspond to stable "stars" with non- 
vanishing scalar field. The final state of the instability 
in this case is an interesting topic for numerical simula- 
tions (see also [l|| and references therein). Note that 
no-hair theorems for positive £ ( in p articular £ > 1/2) 
have been proven yea rs ago [171 . 1 1 81 ] and supported by 
numerical searches [19( . This may imply that either stars 
with £ > 1/6 evolve to a black hole solution in pure GR, 
or that they shed enough mass to leave the forbidden 
region and become a star in pure GR again. 



III. EXTENSIONS 

The LMV mechanism, like other strong-field effects in 
scalar-tensor theories [!, [2(| , relies on a specific coupling 



of the scalar field with matter. We found that an instabil- 
ity is not present for "minimally-coupled" Weyl fermions 
or Maxwell fields, independently of the EOqj. 

Wc have not explicitly shown that the equilibrium con- 
figurations reported here arise as a result of nonlinear 
time evolutions of GR solutions with a "seed" scalar field; 
however, strong indications that this should be the case 
come from numerical studies by Novak in the context of 
spontaneous scalarization Q. It would be interesting to 
perform similar simulations for £ > 1/6. Other possi- 
ble extensions of our work concern the investigation of 
slowly and rapidly rotating stellar models and of their 
oscillation frequencies (see e.g. [HI)- 



IV. CONCLUSIONS 

Wc have reconsidered a generic class of theories where 
a scalar field is nonminimally coupled to the Ricci scalar, 
that were recently shown to give rise to a semiclassi- 



2 The electromagnetic case can be treated as described in 
Ref. This leads to an equation of the form d 2 ^/dr 2 + 

(uj 2 — V) \f = for both polarization states, where we de- 
fined dr/dr, = v7(l - 2m/r), and V = fl(l + l)/r 2 . Mass- 
less neutrino fields yield a similar equation with potential V = 

(2r 2 )- 1 [2kf(r) (k + s/1 - 2m/r) -kry/1 - 2m/rf], where k is 

a separation constant [22J • The potential for Maxwell perturba- 
tions is positive-definite, and therefore no instabilities can arise. 
For neutrinos we were not able to prove stability in general, but 
we did verify that the two specific models for the EOS discussed 
in this paper lead to positive-definite potentials. 
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cal instability. We have pointed out an interesting re- 
lation between the semiclassical instability and the spon- 
taneous scalarization effect in classical scalar-tensor the- 
ories. For certain values of the coupling parameter the 
scalar field can leave observable imprints on the equi- 
librium properties of relativistic stars. Our main find- 
ing is that the LMV instability may provide a "natural" 
seed mechanism to produce spontaneous scalarization, 
reinforcing the relevance of previous studies of compact 
stars in scalar-tensor or f(R) theories of gravitation (see 

e.g. asm). 

We stress that corrections to GR due to scalar fields 
are a generic feature of a large class of unification theo- 
ries. Our work suggests that strong-field modifications to 
GR compatible with weak-field tests may be astrophysi- 
cally viable, with potentially observable consequences in 
the structure of compact stars. It will be important to 
explore the implications of vacuum amplification mecha- 



nisms for tests of strong-field gravity in compact objects 
(see e.g. [26| and references therein). 
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